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ABSTRACT 



o 



Context. The far-IR/sub-mm spectral mapping facility provided by the Herschel-PACS and HIFI instruments has made it possible to obtain, for 
the first time, images of H2O emission with a spatial resolution comparable to ground based mm/sub-mm observations. 

Aims. In the framework of the Water in Star-forming regions with Herschel (WISH) key program, maps in water lines of several outflows from 
yS^ young stars are being obtained, to study the water production in shocks and its role in the outflow cooling. This paper reports the first results of 
C^ this program, presenting a PACS map of the 0-H2O 179 yum transition obtained toward the young outflow LI 157. 

JL Mettiods. The 179 /im map is compared with those of other important shock tracers, and with previous single-pointing ISO, SWAS, and Odin 
i water observations of the same source that allow us to constrain the H2O abundance and total cooling. 

Results. Strong H2O peaks are localized on both shocked emission knots and the central source position. The H2O 179 /jm emission is spatially 

correlated with emission from H2 rotational lines, excited in shocks leading to a significant enhancement of the water abundance. Water emission 

peaks along the outflow also correlate with peaks of other shock-produced molecular species, such as SiO and NH3. A strong H20peak is also 

, observed at the location of the proto-star, where none of the other molecules have significant emission. The absolute 179 yum intensity and its 

7-H ■ intensity ratio to the H2O 557 GHz line previously observed with Odin/SWAS indicate that the water emission originates in warm compact 

^ clumps, spatially unresolved by PACS, having a H2O abundance of the order of lO"*. This testifies that the clumps have been heated for a time 

f^ — long enough to allow the conversion of almost all the available gas-phase oxygen into water The total H2O cooling is ~ 10"' Lq, about 40% of 

7-H the cooling due to H2 and 23% of the total energy released in shocks along the LI 157 outflow. 
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1. Introduction abundance and kinematics in the cold outflow gas components 

(e.g., Franklin et al. 2008, Bjerkeli et al. 2009). Afl these faciH- 
Among the main coolants in molecular shocks, water is the ^^^^^ however, had poor spatial resolution (i.e. greater than 80"), 
tracer most sensitive to physical variations and the temporal evo- ^^^-^^^ ^-^ ^^^ ^jj^^ ^^^ ^^ j^^^^^ ^j^^ ^^-^^ ^f ^1^^ ^^^^^ ^^■^_ 
lution of protostellar outflows, thus representing a very power- ^j^^ ^^^ ^j^^y variations in abundances and excitation within 
ful probe of their shock conditions and thermal history (e.g., individual flows 
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bi Bergin et al. 1998). Water emission and excitation in shocks t 1 • .c 1 1 j^ -im^ a n 1 

Cu ° ... ■ , r , ^ ■ ■ , , T r 1 In this framework, a sample or YSO outflows will be sur 

■ ■ ■ ' were studied extensively tor the hrst time with the Infrared ^ 



Space Observatory (ISO), the first space facility with spectro- 



veyed in different water lines by the PACS and HIFI instru- 

,.,..., ., J J- rr. rn^ j i ments onboard the i/er.sc/!e/ Satellite, as paTt of the key program 

scopic capabilities in the mid- and tar-IR. ISO surveyed the wa- ,,,Tnii /^^r t c^ r ■ ■ ■. rr > irr ^1 ■ 

^ . . . , , r a r 11 WISH (Water In Star-torming-regions with Herscneri). This pa- 

ter emission in a large sample ot outflows from young stellar ,„ ifji^ .• •■ 

,. /,irn^ X -J- 111 • • 1 • jTi • per presents the first results obtained from this survey, consisting 

obiects (YSOs), providing a global statistical picture ot the im- ^„ >,a^o r.u it/-vT i nn i- • .u f 

■' J- ■ , r, 1- J i7 • • ■ ■ of a PACS map of the H2O212-I01 179 yum line covering the out- 

portance ot water in the outflow cooling and ot variations in Its „ „ , ^ rii^^ i • j j • , tt , , 

, J 111 • J / XT- — ,r,r,T flow ot the protostar LI 157-mm, obtained during the //ericneZ 

abundance with shock properties and ages (see e.g., Nisini 2003, „ . ,^ • m t-i nrv i- • i 

T^- 1 1 ^nnA\ f- 11 ■ rn^ 1 ni^rAn j ^j- .c Scicncc Dcmonstration Phase. The 179 um line is the transi- 

van Dishoeck 2004). Foflowing ISO, the SWAS and Odin fa- . .• .u i . 1.11 1 1 r .u n /~v 

.... , . ... , , 1 .. ^jT J 1 tion connecting the lower two back-bone levels ot ortho-HoO. 

cilities made it possible to observe the ortho-HiOtundamental , ■ , s- r , , ■ , ,■ j-i 

,. ^^„ ^,, J- • 1 It IS therefore one ot the brightest water lines expected in col- 

Iine at 557 GHz, providing important constraints on the water ... ,, ., ... ", . •,,, 

^ o ±- Iisionally excited conditions, thus representing an ideal tracer 

. „ ~~~. TT" , .... of the water distribution in shocked regions. LI 157 is a well 

Herschel is an ESA space observatory with science instruments 

provided by European-led Principal Investigator consortia and with im- 

portant partecipation from NASA ' http://www.strw.leidenuniv.nl/WISH/. 
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known outflow driven by a low mass class object (L1157- 
mm, Lhoi ~ 8.3 Lq, D=440 pc, Froebrich 2005). It is considered 
to be the prototype of chemically active flows, given the large 
number of diff'erent species detected in its shocked regions (e.g., 
Bachiller & Perez-Gutierrez 1997). This paper is presenting the 
first of several observations planned for this source by the WISH 
team. 



2. Observations 

Observations were performed on 26 October 2009 with the 
PACS instrument (Poglitsch et al. 2010) onboard the Herschel 
Space Observatory (Pilbratt et al. 2010) in line spectroscopic 
mode, with the grating centred on the H2O 2i2-loi line at 
179.527 /zm. The LI 157 outflow region (of about 6' x 2') was 
covered by 3 individual PACS raster maps, arranged along the 
outflow axis. Each map consists of 3x3 PACS frames acquired 
in steps of 40". The instrument is a 5 x 5 pixel array providing 
a spatial sampling of 9.4"/pixel, while the spectral resolution at 
179 jum is R ~1500 (i.e., ~210 kms"'). The data were reduced 
with HIPE 2.0. Additional IDL routines were developed to con- 
struct a final integrated and continuum-subtracted line map. Flux 
calibrations used calibration files obtained by ground tests that 
remain very uncertain at the time of paper writing, especially 
for extended sources. To evaluate the flux uncertainty, we com- 
pared with the three measurements performed by the ISO satel- 
lite along the outflow (Giannini et al. 2001). To do that, we per- 
formed aperture photometry of the line emission in the PACS 
map within the 80" ISO circular beam. The ratio of PACS to 
ISO fluxes ranges between I.l and 1.8 at the three positions: 
we adopt this as the uncertainty in our quantitative analysis. 
The typical rms noise across the map is of the order of 2 10 *" 
ergs"' cm^^sr"'. 



3. Results and comparison with) othier tracers 

Figure 1 presents the PACS map of the 179 /im line emission. In 
the same figure, the H2O map is overlaid with contours of the 
emission from the H2 0-0 S(l) (Neufeld et al. 2009), CO 2-1 
and SiO 3-2 (Bachiller et al. 2001) transitions. The water map 
exhibits several emission peaks corresponding to the positions 
of previously-known shocked knots, labelled as B0-B1-B2 for 
the south east blue-shifted lobe, and R0-R-R2 for the north west 
red-shifted lobe, following the nomenclature of Bachiller et al. 
(2001). These emission knots represent the actual working sur- 
faces of a precessing and pulsed jet and are thus associated with 
the present location of the active shock regions. With respect to 
CO, H2O emission appears more localized, having a less promi- 
nent diffuse component. About 60% of the total 179 jum flux is 
found within 30" apertures centered on the knots. This could be 
partly related to the line excitation: the 179 /im line excitation 
temperature is ~80 K above the 0-H2O ground state (compared 
to the 17 K for CO 2-1), and the critical density of its upper 
level is above 10*^ cm"^ for T % 500 K. It may however also 
be a consequence of the specific conditions needed to ensure a 
significant production of water The H2O abundance is indeed 
significantly higher only in shocks strong enough to release the 
water ice located on grain mantles by sputtering and grain-grain 
collisions or to activate the gas-phase reactions that convert the 
gas-phase oxygen into water. Both these processes become ef- 
ficient at shock velocities v, k, 15 kms"' (Caselli et al. 1997, 
Jimenez-Serra et al. 2008, Kaufman & Neufeld 1996). In this 
respect, we note that the H2O emission peaks correspond rather 



closely to both the position and the relative intensity of the H2 
rotational emission (with the H2O 179 yum/H2 17 //m ratio in the 
range ~(2-3)xlO"^foralltheH2 peaks). Peaks of low-/ H2 pure 
rotational lines are associated with warm gas (with T ~ 300-500 
K) excited in low velocity non-dissociative shocks that are trac- 
ers of regions in which a high H2O abundance is expected. Other 
molecules are known to have strongly enhanced abundances in 
shocks. One of the most well studied of these molecules is SiO, 
for which Fig. 1 shows that, like water, its emission is very lo- 
calized around the shocked knots. A similar behavior is found 
for other molecules, such as NH3 and CH3OH (Bachiller et al. 
2001, Tafalla & Bachiller 1995). 

The strongest water peak is located at the position of the B 1 
knot, which is known to be the most chemically active of the 
LI 157 spots (e.g. Bachiller & Perez Guitierrez 1997, Benedettini 
et al. 2007, Codella et al. 2010). This knot at near-IR wave- 
lengths appears as a bow shock with intense H2 2. 12//m emission 
(Davis & Eisloffel 1995) and has a significant H2 column den- 
sity enhancement (Nisini et al. 2010). Although the spatial res- 
olution of the present observations prevents us from completely 
resolving the bow shock structure, the observed morphology at 
the Bl/BO positions suggests that water emission is mainly lo- 
calized at the bow apex and eastern wing. A similar morphol- 
ogy has been observed for molecules such as SiO, NH3, and 
CS (Benedettini et al. 2007, Tafalla & Bachiller 1997), while 
other shock produced molecules, such as CH3OH, noticeably 
have emission localized on the bow western wing (e.g. Codella 
et al. 2009). This behavior probably relates to an asymmetry in 
the excitation conditions along the bow structure, most likely in- 
duced by the jet precession or the propagation of shocks in an 
inhomogeneous medium. 

Strong, spatially unresolved, water emission is also detected 
on-source. This localized emission can originate in different 
components, including shocks impacting on a dense medium at 
the jet base, the infalling protostellar envelope, or emission from 
a UV-heated outflow cavity, as discussed in van Kempen et al. 
(2010) for the HH46-IRS case. The precise origin of this emis- 
sion will be investigated by dedicated Herschel observations, but 
we note here the interesting evidence that no other molecule ex- 
hibits significant emission at the central position. In particular, 
the non-detection of strong emission from molecules such as 
CH3OH indicates that grain ice mantle evaporation in the proto- 
stellar envelope is unlikely to be the origin of the on-source H2O 
emission, since the two molecules should desorb at similar tem- 
peratures. The non-detection of the H2 0-0 S(l) line at the cen- 
tral position is also remarkable. This may be caused by the heavy 
extinction close to the central source. Assuming an intrinsic 
H2O 179/im/H2 17yum ratio in the range of that observed along 
the outflow, we estimate that A,, on-source should be k, 150 mag 
to be able to explain the H2 line non-detection. Alternatively, C- 
type shocks with very high pre-shock densities (>10^cm"^) and 
velocities between 20 and 40 km s"' are expected to have a large 
H2O/H2 cooling ratio (Kaufman & Neufeld 1996). 

4. Water abundance and totai cooling 

To constrain the range of water column densities that could pro- 
duce the observed 179 /im emission, we consider the SWAS and 
Odin observations of the H2O lio-loi 557 GHz (538/im) line 
observed in this outflow (Franklin et al. 2008, Bjerkeli et al. 
2009). Given the large size of the apertures of these two instru- 
ments relative to the PACS spatial resolution, we evaluate here 
only properties averaged over large outflow regions. In partic- 
ular, we consider the Odin observations acquired towards the 
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Fig. 1. Continuum subtracted PACS map of the integrated H2O 179 //m emission along the LI 157 outflow. Offsets are with respect 
to the L1157-mm source, at coordinates q'(2000) = 20:39:06.2, 5(2000) = +68:02:16. The different emission peaks are labelled 
following the nomenclature adopted by Bachiller et al. (2001) for individual CO peaks. The same map is shown in the other panels 
with overlays of other tracers, namely H2 0-0 S(l) at 17;um (Neufeld et al. 2009), CO 2-1, and SiO 3-2 (Bachiller et al. 2001). The 
spatial resolution of these images are ~H", for H2 and CO, and 18" for SiO. Note that the H2 observed region does not cover the 
B2 and R2 shocked peaks. 



blue (B) and red (R) outflow lobes at offsets (-h29",-52") and 
(-21",+ 121") (BjerkeH et al. 2009). The 179A/m/557GHz in- 
tensity ratios are obtained by diluting the PACS observations 
to the 126" Odin resolution. The same procedure was adopted 
for the SWAS observation that encompasses almost the entire 
LI 157 PACS mapped region with its 3.5'x 5.0' elliptical aper- 
ture. Figure 2 presents large velocity gradient (LVG) predictions, 
assuming a slab geometry, of the 179 jum line brightness ver- 
sus the 179//m/557 GHz line ratio, compared to the observations 
combined above. The absolute brightnesses are those averaged 
within an area enclosing 90% of the total PACS emission inside 
each considered Odin/SWAS aperture. These emitting areas are 
5.9x10 ^ 8.0x10 ^ and 2.7x10"^ sr for the R, B, and the SWAS 
apertures, respectively. The line intensity derived in this way was 
considered to be a lower limit to the true 179 //m brightness if 
the PACS emission originates in a clumpy medium, of which 
the clump size is smaller than the Herschel diffraction limit at 
179 jum. 

In the figure, observations are indicated as boxes that take 
into consideration the uncertainty of a factor of about 1 .5 in the 
179 /vm flux, estimated by comparing with the ISO observations 
(Sect. 2). Theoretical curves were derived as a function of the o- 
H2O column density, using the RADEX code (Van der Tak et al. 
2007) assuming temperature and density conditions measured 
from the H2 Spitzer observations or ground-based millimeter 
observations (Nisini et al. 2010, Nisini et al. 2007, Mikami et 
al. 1992). The temperature is between 300 and 500 K and the 
density is in the range l-5xl0^cm"^, the blue lobe being on av- 
erage colder and denser than the red lobe. Part of the 557 GHz 
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Fig. 2. LVG theoretical predictions of the 179 yum line brightness 
versus the 179/im/557GHz line ratio, compared with observed 
values. See text for the details. 



emission can arise from a gas colder than these assumed values, 
given the lower excitation temperature of this line with respect to 
the 179 yum line. To evaluate the effect of different temperature 
components along the line of sight on the ratio of the two con- 
sidered transitions. Fig. 2 also plots the theoretical predictions 
assuming a temperature stratification where the column density 
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in each layer at a given T varies as T * (Neufeld & Yuan 2008). 
A minimum and maximum temperature of 100 K and 4000 K, 
respectively are assumed, and b values between 2 and 4, i.e., 
the range of values that consistently fit the H2 rotational lines 
(Neufeld et al. 2009). These curves give the same range of pre- 
dicted values as the single T curves, indicating that contributions 
from high-temperature gas do not significantly affect the consid- 
ered transitions. 

Several general conclusions can be drawn from the inspec- 
tion of Fig. 2. Firstly, the data are consistent with model predic- 
tions only if we assume that the real emitting areas are smaller 
than those estimated from the PACS map. In particular, agree- 
ment with the theoretical curves is found for covering factors (/c) 
~ 0.1-0.2, which suggests that the emission is concentrated on 
some unresolved emission knots that together do not fill an area 
larger than a few tens of arcsec. This is not unexpected, since 
interferometric mm observations illustrate the extreme dumpi- 
ness of the shocked gas, individual knots being of sizes of a few 
arcsec each (e.g., Benedettini et al. 2007, Lefloch et al. 2010). 
We note that the typical length scale for planar C-type shocks 
at the considered densities is of the order of lO"" cm, i.e., about 
1/10 of the PACS spatial resolution at D=440 pc. The observed 
179//m/557GHz ratios, ranging between 10 and 20, are con- 
sistent with A^(H20) ~ 2-9x10^^ cm"^ (assuming a Av = 15 
km s ' from the 557 GHz line width). The H2 column densities, 
averaged within the PACS emitting areas, were measured from 
the H2 mid-IR rotational lines and results in ~ 5xlO'^cm"^ in 
both regions covered by the B and R observations. The water 
abundance in the unresolved clumps is therefore estimated to be 
~ A^(H20)/A^(H2)x/, -0.6-3 lO""* (with a H2O o/p ratio of 3). 
Table 1 reports in more detail the range of values derived in each 
considered aperture. The total mass of the shocked gas involved 
in the 179 yum emission is of the order of 5 lO"-' Mq, which 
is only a small fraction (~ 1/100) of the total mass of the out- 
flow estimated from CO observations (e.g. Bachiller et al. 2001). 
Lefloch et al. (2010) show that H2O components with different 
velocities are discernible in the 557 GHz data acquired by HIFI 
in a 40"beam centred on the L1157-B1 knot. They separately 
analyse the different velocity components, confirming that small 
filling factors are required to explain their observations and find- 
ing that the component of higher velocity is the one exhibiting 
the water abundance of the order of 10 '*. Lower H2O abundance 
values, between lO"*" and 10"^, were estimated using only the 
SWAS and Odin 557 GHz emission, assuming that the 557 GHz 
emission originates in the same cool gas traced by the low-/ CO 
emission, thus a gas with a larger covering factor and lower tem- 
perature than considered here (Neufeld et al. 2000, Franklin et 
al. 2008, Bjerkeli et al. 2009). Combining ISO-179yum emission 
and SWAS observations, Benedettini et al. (2002) derived a wa- 
ter abundance for the warm shocked gas of ~ 5 10"^, thus in the 
lower range of values estimated in the present analysis. However, 
the ISO observations did not cover the entire LI 157 outflow 
179 yum emission, and the inferred ISO 179yum/SWAS 557 GHz 
ratio was underestimated by about a factor of 2. Given the con- 
sidered conditions, the 179 //m line contributes to about 30-40% 
of the water emission in the outflow: the total estimated H2O lu- 
minosity is ~ 8-9xlO"^L0, which is about 40% of the total H2 
shock luminosity (0.2 Lq, Nisini et al. 2010) and about 23% of 
the total shock cooling in the LI 157 outflow, if we also consider 
the contributions given by CO and [O i] derived from ISO ob- 
servations by Giannini et al. (2001). The high water abundance 
estimated in the present analysis is consistent with predictions 
of non-dissociative shock models, in which water is mainly pro- 
duced by endothermic reactions, activated at T ^ 300 K, where 



Table 1. Estimated water abundances 



K 



"(Hz)" 



lV(oH^O) XCHaO)* W 

10'*cm-2 10-* sr 



B 
R 



14-20 
10-14 



300 
500 



TkF" 

1 10^ 



3-9 

2-4 



0.6-3 
0.8-2 



■ 5.2 10-" 

■ 9.4 10-' 



" 179yum/557 GHz ratio witliin the Odin aperture. 

* See text for references and assumptions on T, n and A'(H2). 

^ Effective emission area that reconciles the observed and predicted 

179 yum line intensity within the Odin aperture. 



all the available gas-phase oxygen is converted into H2O, or by 
the sputtering of icy grain mantles behind the shock. According 
to Bergin et al. (1998), the time needed to complete this process 
is of the order of 10^ yr, for T - 400 K. This is comparable to 
the shock timescales estimated from H2 observations of individ- 
ual emission knots of the LI 157 outflow (Nisini et al. 2010), thus 
supporting the idea that the water in this outflow has had time to 
reach its maximum allowed abundance. 



5. Conclusions 

We have presented a PACS spectral map of the H2O 179 /im 
transition obtained toward the LI 157 protostellar outflow. Strong 
water emission peaks have been found at the location of 
previously-known shocked spots and correlate well with H2 mid- 
IR rotational lines, as well as other important shock tracers, such 
as SiO and NH3. The absolute 179 yum intensity and the in- 
tensity ratios with respect to the previously-observed 557 GHz 
line, indicate that the water emission originates in warm compact 
clumps, spatially unresolved by PACS, that have a H2O abun- 
dance of the order of 10"'*. The total H2O cooling has been esti- 
mated to be of the order of 8-9 lO-^L©, representing about 40% 
of the cooling due to H2 and 23% of the total energy released in 
shocks along the LI 157 outflow. 

Additional Herschel PACS/HIFI observations of the LI 157 
outflow are planned by the WISH program. These will enable 
us to investigate variations in the water abundance within the 
outflow and correlate these with kinematical information. 

Acknowledgements. This program is made possible thanks to the HIFI guaran- 
teed time and the PACS instrument builders. 
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